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Bose-Einstein condensation of a relativistic ideal Bose gas in a rectangular cavity is studied. 
Finite size corrections to the critical temperature are obtained by the heat kernel method. Using 
zeta-function regularization of one-loop effective potential, lower dimensional critical temperatures 
are calculated. In the presence of strong anisotropy, the condensation is shown to occur in multisteps. 
The criteria of this behavior is that critical temperatures corresponding to lower dimensional systems 
are smaller than the three dimensional critical temperature. 



I. INTRODUCTION 



Bose-Einstein condensation (BEC) first predicted in 1924 |ï]-[i| has become an explosive field of research in recent 
years Due to the temperature dependence of the chemical potential in a Bose-Einstein distribution function, 

macroscopic number of bosons start accumulating onto the ground state at the critical temperature. This phenomenon 
is typically phrased as the phase density of partides being large enough such that all partides characterized by the 
de Broglie wavelength overlap to form a condensate. Liquid helium, which becomes a supernuid at the transition, 
has until recently been known as the only substance which shows this behavior. However, strong interaction between 
V/ helium àtoms has been an obstacle for a complete understanding of the mechanism of condensation. 

Recent technological progress in atom cooling techniques made it possible to achieve Bose-Einstein condensation 
for neutral àtoms. Weakly interacting nature of these àtoms enables one to understand the condensation within the 
perturbative framework. It also provides an ideal testing ground for some fundamental aspects of quantum mechanics 
in a controlled environment. Studies of similar phenomena, stimulated by ràpid progress in this subject, are no 
longer restricted to condensed-matter physics and atom-optics but start to involve other areas in physics, such as 
nuclear-particle and astrophysics ||. 

In this paper we study the effect of a finite size container on the condensation. For a finite size system, the absence of 
thermodynamic limit alters various critical behaviors defined and expected for a bulk system |)|-|lïj . Thermodynamic 
Oh, quantities such as the free energy has a surface term which vanishes in the thermodynamic limit, causing a shift in 
' the critical point. Finite size effects in Bose-Einstein condensation manifest themselves as the rounding-off of the kink 
in the specific heat at the critical temperature^. Off-diagonal long range order stays in a finite range at the critical 
point. 

q-( For a finite size quantum system, the invariant operator of small fluctuations has a discrete spectrum. The relevant 
• • i dimensionless parameter r/i which characterizes finite size effects near the transition is given by r/i = XgdB / Li for a 
nonrclativistic system, where Xbúb is the thermal de Broglie wavelength and Li(i — 1,2,3) are the system sizes in 
the three spatial directions. r/i — (3 /Li for a relativistic system, and rji — (3lüí for a harmònic oscillator with natural 
$H ' freqüències u>i(i = 1, 2, 3).f] 

Presence of anisotropy adds more variety to the critical behavior. Suppose rji > 1 for some i (for a nonrelativistic 
case, this implies Li < Xeds) then only the lowest mode in the i-th direction contributes significantly to the dynamics 
of the system. The motion in the i-th direction is frozen out and the system has an infrared behavior effectively 
equivalent to a system with one less dimension in that direction |Ï2-17|. Thus we can classify the dynamics with an 



effective infrared dimension (EIRD) into the following four cases dependent on the degree of anisotropy: 
Case 1; íyi , í?2 , ?73 > 1 — * EIRD = 0, 
Case 2; 771, 772 > 1 > rj 3 — > EIRD = 1, 
Case 3; 771 > 1 > 772,773 — > EIRD = 2, 



*E-mail address: kshiok@phys.ualberta.ca. Present address: Theoretical Physics Institute, University of Alberta, Edmonton, 
Alberta T6G 2J1, Canadà 
^E-mail address: hub@physics.umd.edu. 

1 The bulk specific heat of an ideal Bose gas is not divergent, but shows the discontinuity in its derivative. 

2 In this paper, we use the units where fes = h = 1, which renders rji dimensionless. The results in ordinary units can be 
obtained by replacing u) here by hui and T here by fcsT. 
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Case 4; 1 > 771, 7j 2 , V3 — > EIRD = 3. 
In this paper, we mainly study Case 4 where modes in all three directions are excitable. As the temperature is 
lowered, the crossover behavior between higher- and lower-dimensional excitations can be observed. Each mode is 
labeled by three quantum numbers associated with the excitation energy in each direction. In the presence of strong 
anisotropy, with these quantum numbers, it is meaningful to split the whole particle spectum into either zero, one, 
two, or three-dimensional excitations. The ground state, being the state with the lowest quantum numbers, is viewed 
as a zero-dimensional excitation. Let us denote the number of modes excited in the corresponding directions as 
N , Nx, N 2 , N 3 , respectively. 

We can define an n-dimensional critical temperature T n £> (n = 1,2,3) as the temperature at which all the n- 
dimensionally excited modes are saturated: 



where the thermodynamic limit in each case is taken differently. As we set the total number of partides N — > 00, 
we obtain T 3 u by tuning rjx, i] 2 , and 773 to zero, T 2 n by tuning 772 and 773 to zero while 771 fixed, and T\u by tuning 
773 to zero while r\\ and 772 fixed. Finite size corrections necessarily modify the above definitions, since they involve 
excitations in lower dimensions. In Section 3, we will discuss this aspect in detail. 

By changing the edge lengths of a cavity or oscillator freqüències for a magnètic trap, it is possible to control the 
critical temperature and realize the lower dimensional condensation. In particular, when T\r> < T 2 £> < T 3 o holds, 
condensation is expected to oceur in three steps: As the temperature is lowered, condensation into two-dimensionally 
excited modes begins at T 3 d when three-dimensionally excited states saturate. At the two dimensional critical 
temperature T 2 £> condensation into one-dimensionally excited modes begins. The condensation onto the ground state 
does not oceur until one dimensional critical temperature Txd is reached. 

In a finite size system, the reduced chemical potential e = (3(Eq — /i) does not vanish. From the expression for the 
ground state contribution, we can still assume e ~ up to order 1/Nq. This condition is justified as long as Nq is 
close to the total number of partides, or equivalently, the temperature is lower than the critical temperature. 

Although work on BEC in rclativistic systems has a long history, modern treatment using quantum field theory did 
not begin until 1980's |p^| , ^9| , pT| -[23]| . At relativistic temperatures, T > 2m, where m is the mass of the relativistic field, 
pair creation-annihilation effects become nonnegligible, and the particle number is no longer conserved. However, the 
global U(l) gauge symmetry of the Hamiltonian guarantees the existence of a conserved charge based on Noether's 
theorem. The net charge Q in relativistic field theory is given by 



Partides and anti-particles have chemical potentials opposite in sign due to the fact that they carry opposite charges. 
Taking this fact into account, the positive-definiteness of the particle number of the partides and antiparticles with 
energy En requires that |/x| < m. 

Another important point is the relation between spontaneous symmetry breaking (SSB) and Bose-Einstein conden- 
sation jÏ9| . Condensation into the ground state results in a nonzero vacuum expectation value of the field. Hence 
BEC can be interpreted as a SSB of the local gauge symmetry. This argument presumes that the chemical potential 
reaches its critical value at the critical temperature. However, for a finite system, this is generally not the case. 

Bosc-Einstcin condensation of a relativistic noninteracting quantum field in a rectangular cavity is studied in this 
paper. Similar aspeets for a nonrelativistic ideal Bose gas in anisotropic magnètic traps is discussed in a companion 
paper pp[ |. In Section 2, we derive the effective action which includes one-loop quantum corrections to the classical 
action and use the £-function regularization to evaluate the grand canonical thermodynamic potential p4]~p7|. The 
generalized £-function is written in terms of 0-functions via a Mellin transformation. We then use the asymptotic 
expansion of 0-functions for a large system size and a small massive field to see the finite size correction to the total 
charge and the critical temperature. This asymptotic expansion is a special case of the more general class of short 
time expansion of the heat kernel which is used for spectral analysis on an arbitrary differentiable manifold ps} ] . The 
terms in the expansions correspond to the volume (Weyl), area, and edge contributions, etc |2^ , ^9| , ^o| . In Section 3, 
we consider the effects of accidental degeneracy in a discrete spectrum and show that the highly oscillating behavior 
of the density of states is large enough to dominate over the higher order terms in the asymptotic expansion. We 
introduce an infrared cutoff to include the lowest mode contribution properly and estimate the lower-dimensional 
critical temperature accurately. In the last part of this paper, we discuss the multistep behavior of condensation 



3-dimcnsional temperature; N 
2-dimcnsional temperature; N 
1-dimensional temperature; N 



N 3 (T 3D ), 

N 3 (T 2D )+N 2 (T 2D ), 

N 3 (T 1D ) + N 2 (T 1D ) + Nx{T lD ), 



(1) 
(2) 
(3) 




(4) 
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process in the presence of strong anisotropy pjmlfl . The conditions for one, two, and three-dimensional condensations 
are clarified. The relevant critical temperatures are obtained. 

As discussed and observed in (3|^^1 for a weakly interacting gas, the corrections to bulk ideal-gas ground state 
occupation number and critical temperature are well-explained by the finite size effects. Interaction effects on thosc 
quantities are negligibly small. Hence we expect the results discussed in this paper will still hold for weakly interacting 
gases. On the other hand, interaction effects are known to afïect higher moments such as the specific heat significantly 
and considered to be essential in explaining the observed specific heat data. Extension of our analysis to the strongly 
interacting case is a nontrivial problem which deserves further careful study. 



Bose-Einstein condensation of a relativistic gas could be relevant to cosmology in the dark matter problem |34|, or 
for inflationary universe pq| . Our problem is directly related to condensation of positronium in a cavity discussed 
in p6| . Although we restrict our study to a rectangular cavity similar results are expected for systems with a finite 
boundary. In the presence of an irregular boundary the dynamics of the system can be nonintegrable or chaotic. The 
implication of our result for nonintegrable systems of this kind is of particular interest [ |37| . The extension to curved 
spacetimes can also be obtained by similar mcthods pq HW. 



II. EFFECTIVE ACTION AND HEAT KERNEL 
A. One Loop Effective Action 

The action of a free complex scalar field 

S Q [cf\ = J éx^^d^ - m V0] (5) 

is invariant under the global U(l) gauge transformation 

(f>(x) -> e"V(x). (6) 

The corresponding Noether current is 

J^x) = itfd^ - id^cp, (7) 

with the total charge 

Q= í d 3 x,J (x)=i i d 3 x[tf(j) (8) 
Jv Jv 

The integration is over the volume V of the cavity, here assumed to be rectangular of edge lengths Li(i = 1,2,3). 
Decomposing </>(x) into real and imaginary parts such that cf>(x) = ï/ \/2[<f>i(x) + i<f>2{x)], Eq. (|5|) becomes 

SqW\ = \J dïx^^Mi + 8*48^2 - m 2 4fa - m 2 ^}. (9) 

The Hamiltonian for this action is 

H= \ I d^ln 2 + tt 2 2 + + (V0 2 ) 2 + m 2 ^ + m 2 4fo], (10) 

^ Jv 

where ~k\ — 4>i,tt2 — $2 are the momentum fields canonically conjugate to (f>i and <j>2 ■ The total charge becomes 

Q = J d 3 x(4>2Tïl - 0l7T 2 ). 

The grand canonical partition function for this system when brought in contact with a heat bath at temperature 
T = 1/(3 is given by 

Z = Tre- f, &-"à\ (11) 

where H and Q are the Hamiltonian and the total charge operators respectively and \i is the chemical potential. Eq. 
( |ï"ï| ) in Hamiltonian form has a path integral representation 

Z = [ DnD<l>exp[ [ drí d 3 x[in 1 4 + in 2 (f> 2 - H + /í(0 2 tti - 0itt 2 )]], (12) 
J Jo Jv 
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where fa = d T fa. In the spirit of (imaginary time) finite temperature field theory, a periòdic boundary condition is 
imposed on fa, with fa(0, x) = fa(P,x). We perform an integral over the momentum field and obtain 

Z = J D^e'^, (13) 

with the action 

S[fa = dr d 3 x[^(fa - iiifaf + l -{fa + líifa) 2 

+ \(yfa? + \{Vfaf + \m\fa 2 + fal)]- (14) 

Using the background field decomposition cf> — fa + ip with fluctuation <p, and expanding the action in Eq. ( |Ï4| ) 
around the classical solution fa which minimizcs the action 



S[fa = S[fa] + \Y. aSt^ + °^ 3 )- ( 15 ) 



The partition function can be written as 



where T[fa] is the effective action and Ay[0 c ] = 5 2 S[fa,}/Sfaí 
The functional measure Dip is defined as 



where c n are the coeficients of an eigenfunction expansion of ip and / is a constant with unit of length. Then the 
functional integral in Eq. (^) can be evaluated as 

i r°° 

TT-=-/ dc n e-^^ ^mt^K^fa])- 1 ' 2 . (18) 

„ V27T/ J_oo 



The effective action to one loop order is given by [421 



Y[fa] = S[fa] + ~ logDet(l 2 A«[^ c ]). (19) 



The second term in Eq. ( |Tgj ) can be split into two parts as 



logDet(; 2 A y [^ c ]) = logDet(í 2 A + [(/. c ]) +logDet(; 2 A_[^ c ]), (20) 



where A± = — (<9 T ± fï) ~ V + m . The eigenvalues of A± are given by 



where ujn is the eigenvalue of — V 2 



B. £-function Regularization 

The generalized £-function for an elliptic differential operator O is defined by 

Co( S )-TrO- s -^A w s , (22) 



N 
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where Atv are eigenvalues of O. From Eq. Q22J), 

CÓW = -EVlogV- (23) 

JV 

Thus 

logDet(/ 2 ü) = Co(0)logZ 2 - Có(0). (24) 
Using a Mellin trasformation defined by 

A- s = — / dr r*"^-^ (25) 

we can write the generalized £-function for Ajj as 

C A± (s)=TrA ± - s =^A ± ^ 
1 í°° 

where A±j\r and K\ ± (t) are eigenvalues and the heat kernels for A±. 
Here 





n= — oo TV 

yielding 

Ca(.s) = Ca+(s) = Ca» 



(26) 



oo _ 

K A± (r) = K (t) = J2 E ex P[-^[(^+^) 2 +^/v + m 2 ]] (27) 



drr-^o^r), (28) 

r(s) 7 

where /3 = /3/2vr |ï]]4| and 

^ (/3 2 r) - A^Me-^-^^^s^Kr/Tr), (29) 
where 9^(z\t) = 1 + 2J2^=i e% ' KT cos(2nz) is a #-function |Ï6j. The heat kernel K(t) for —V 2 is defined by 

K(t) =^2e- 02üJNT . (30) 

TV 

With this, the effective action can be expressed in terms of Ca(s) as 

ry> c ] = ^[^ c ] + CA(o)iogZ 2 ~C/;(o). (3i) 

We first consider Neumann boundary conditions at the boundary of the cavity. The corresponding eigenfunction 
for -V 2 is 



2 sTmiX! ( Tm 2 x 2 , ,Tm 3 x 3 

<M X ) = \j cos(— — ) cos(— — ) cos(— — ) (32) 

L·\L·2L·z L·i L2 L·z 



and the eigenvalue wtv is 



^(^) 2 + (^) 2 + (^) 2 , (33) 

L·i L 2 L·j, 



where m = 0, 1, 2, • • • (i = 1, 2, 3). 
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The eigenfunction for Dirichlet boundary conditions can be written similarly as 



0n(x) = \l LiL2Ls sm( -~^~) Sln( -~L^~' ) Sm( -~^~^ *- ' 



and the eigenvalue u>n is 
where m = 1, 2, • • • (i = 1, 2, 3). 



^v = (^) 2 + (^) 2 + (^) 2 , (35) 

L\ L 2 L 3 



C. Asymptotic Expansion of the Heat Kernel 



The heat kernel for all accessible quantum states is given by 
K{t) = e-P 2uJNT 



N 
oo 



E2 2 2 x -S 2 2 2 x ^ 

e - Vl TT n x T e -7J 2 7T U 2 T g - 

"l ri 2 n 3 

= ^[^(OIVttt) ± lp 3 (0|i % 2 7rr) ± lp 3 (0|i% 2 ^) ± 1], (36) 

where 77^ = /3/ii for i — 1,2,3. Positive (negative) signs correspond to Neumann (Dirichlet) boundary conditions. If 
we assume Li » (3 or r\i << 1 for % = 1,2, 3, we can make use of the asymptotic behavior of 0-function as r — * 

03(O|ir) - 4= ( 3? ) 
to obtain the following asymptotic property for the heat kernel 

K(t) -> h 1 ± 1][ 1 ± 1][ 1 ± 1] 
» ^Jri^TtT V 7 /^ 7 "" V^kt 

1 1.1 1 1 , 1.111,1 

± — ( + + ) + 3-7=(— + — + — ) ± r 



8(7rr) 3 / 2 77i772í73 87rr 771772 772773 7737/1 8v/7rT 771 772 773 
A 3 A 2 A 1 



P3 t 3/2 p2 T p T l/2 



A , (38) 



where 



±A = 1/8, 

Ai = Li + L 2 + L 3 /8tt^ 2 , 
±A 2 = L X L2 + + LaL^Sn, 

A 3 = L^Lz/^ 2 . (39) 

Note that the leading finite size correction A2 has opposite signs for Neumann and Dirichlet boundary conditions. 
This fact results in the opposite shift in the critical temperature as we will see in Eq. (|5^), where the coefhcient 6 3 / 2 
is proportional to A2. This expansion is equivalent to the short time expansion of the heat kernel used in spectral 
analysis on an arbitrary Riemmanian manifold. The first term in Eq. ( |3^ ) is the Weyl term, the second term is the 
boundary contribution, etc p9|j30|] . Since the model is integrable, the heat kernel is factorized into contributions from 
each dimension. 

Using the above expression in Eqs. ( p8| ) and (|2^), we obtain the generalized £ function for A,j 

f52s 3 A roc 

1 v s ; ,_ n P Jo 



k=0 

2s 



vi* r 



dT r s-fe/2-l e -(m 2 -A 1 2 )^ 2 r 



n =l gr=o ^ 



G 



Let us denote the first term in Eq. (^0|) on the right hand side of the equality as Ci( s )j the second term as £2(3). 
Then 



and 



= E r{S ~^ /2) Mrn 2 -^) k/2 - s (41) 



v ' fc=0 H J ° n=l ç=0 v 

(-l)P(m 2 - /i 2 )?/? 2 *^ ^ _ n2r « (-l)9(2n/í/?r) 2<! 



x 



^ p! ^ ^ (2«)! 

p=0 r n=l q=0 K H > 



oo2s 3 A 00 00 
^ > k=0 P p=0 q=0 

x t^(m 2 -M 2 ) P /3 2ï, (2M/3)^. (42) 
We expand Eq. (^), assuming P/Li, mj3 (/i/3), and (to 2 — /i 2 )/3 2 are small quantities 3 , and obtain 

C2(S) = ï|y [ F^ 7/2 ^ 4 - 2s ^ 2 - s ) + ^f- 2s " 5/2 C(3 - 2a)r(3/2 - a) 

+ /3~V 2s - 3 / 2 C(2 - 2a)r(l - a)[Ai - A 3 [m 2 + (2a - 2)^ 2 ]] 
+ C(2s)r( S )[A, - A 2 [m 2 + (2a - l)^ 2 ]] 

+ /3C(2s + l)r(s + l/2) (43) 
x [A 3 [(m 2 - M 2 ) 2 /2 + (a + 3/2) [s + l/2)2 M 4 /3] - A^m 2 + 2s/i 2 )]]. 

The one-loop efïcctivc action then has the form 

r[&] = 5[^]-c;(o) + CA(o)io g ; 2 

+ S + + f + Cl/2 + C ° + c -V2^1og^ + • - , (44) 



where 



* 1/2 A 

c 3 = -—A 3 , 



C2 



C(3)A 



2 



ci = —{A 1 -A 3 [m 2 -2 l· i% 



ci /2 = 2[A -A 2 [m 2 - /i 21 



For large size (Li » j3) and small mass m « T, we do not need to assume higher temperature. See also 
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co = 21og(2^4 - ^(™ 2 - l· i 2 ) 1/2 A 1 - [(to 2 - M 2 )log(TO 2 - M 2 ) 

+ 21og(2^)TO 2 - (21og(27r) + 1)^2 - ^(m 2 - ^ f^A^ 
C-i/2 = -^- 1/2 [A 3 [(to 2 - M 2 ) 2 + /i 4 ] - 2A im 2 }. (45) 
The total charge can likewise be written as 

-i dr[<t>) 



Q = 



(3 du 

hiT 2 + b 3/2 T log T + b l T + b 1/2 log T + b + ■ ■ ■ , (46) 



where 



0-2 = — - — A 3 , 



h/2 = 4/uA 2 , 

h = - M [4^ 1/2 (to 2 - fi 2 ) 1/2 A 3 + 2[2 + 1o S (to 2 - n 2 )]A 2 
2tt^Ai 2 Ap 

(TO 2 — /i 2 ) 1 / 2 TO 2 — /i 2 ' 

27TV2' 

P ln\ , o.A , 32 ^ ..2 q A 1 -,-72n 



2 



bo = 2^72 [2CM1/2) + 3 7 ) + y^ 2 - 8^ - log(Z 2 )C], (47) 

where C = A 3 (4fi 2 — to 2 ), 7 is the Euler constant, and ip is the Digamma function (?/>(l/2) = —7 — 2 log 2 = 
— 1.96351 • • ■). Note that b 2 gives the bulk term discussed in |jïs|| . 

The above derivation is based on an asymptotic expansion of the heat kernel which assumcs a continuum spectrum. 
The same assumption may not be justified in low dimensions where the density of states does not increase as rapidly 
with energy due to the restricted degrees of freedom. Thus, the continuum spectral approximation has to be modified 
accordingly. Indeed this type of expansion does not reproduce the bulk term which appeared in |Ï8|| in the two 
dimensional case p0[ . Also it is not straightforward to define one dimensional critical temperature T\r> with the 
method described above due to the existence of inverse powers of to 2 — /i 2 in Eq. (|47|) . It is known that the chemical 
potential does not reach its critical value for a finite system at the critical temperature. For a relativistic field theory 
in curved spacetime, this aspect has been studied in 

III. FINITE SIZE EFFECTS AND MULTISTEP CONDENSATION 

A. Finite Size Effects and Discrete Spectrum 

In this section, we use an alternative method to treat the discrete spectrum in a more appropriate way. Rewriting 
the heat kernel for all accessible states in Eq. (prf) gives 



E2 2 2 x 2 2 2 x 2 2 2 
g-í/i ir % ry e -> )2 7T n 2 r \ ^ e ~v 3 * «3 r 

m n 2 113 

00 00 00 



ni n 2 n 3 



where qi — e Vi 17 r . Throughout this section, we assume that L 3 is an integral múltiple of L\ and L 2 , such that 
a\L\ = a 2 L 2 = L 3 for some integers a±, a 2 .R K(t) becomes 



Integer assumption here is not essential but for calculational convenience. 
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oo oo oo 

2 2 i 2 2 i 2 

1 a 2 n 2 +n 3 



ni n 2 n 3 
oo 

= X> 3 (n)í 3 », (49) 

where r3 (n) is the number of solutions of the Diophantine equation n = a 2 n 2 + a^n^ + n 3 in natural numbers [ Esl . 

Let us define the function A^e) which counts the number of points with integer coordinates inside the ellipsoid 
whose x,y,z-intercept are s/e/ai, ^flja-i, ^ft, respectively.^ In Appendix.A, we give a derivation of the exact formula 
for Md{£) for arbitrary dimension d which gives the cí-dimensional cumulative density of states A/"d(e) in the way 
described below. 

For d = 1, 2, 3, they are related to each other by: 

Af 3 ( £ ) = 8AÍ 3 (e) ± l2AÍ 2 (e) + 6A?i(e) ± 1, 
AT 2 (e) = 4Ar 2 (e)±4Ar 1 (e) + l, 

M(e) = 2AT 1 (e)±l, (50) 

where the upper (lower) signs correspond to Dirichlet (Neumann) boundary conditions. 

Inverting (|5C|), we obtain the expression for Afd(e) (d = 1, 2, 3) in terms of A/rf(e) (<i = 1, 2, 3) as 

Aí 3 (e) = i [jV 8 (e) T 3AA 2 (e) + 3M(e) T 1] , 

AT 1 ( £ ) = i[M( £ )Tl]. (51) 
From Eqs. ( |5Ï| ) and (A4), we readily obtain 



rr f s n 7T£ r 1 1 1. A , , . . 

A/ 3( £ )=fi T "H" + — +— +A(e). (52) 

d oia2 8 aio,2 ai a 2 

The first term in Eq. (^2|) simply comes from the volume of the ellipsoid, the second term originates from compensating 
the oversubtracted points on the three coordinate planes. The residual term A(e) includes, in addition to the terms 
corresponding to higher order contributions in the asymptotic spectral expansion, the error of approximating cubes 
located on the surface of the sphere by a smooth surface. This error is ascribed to what is known as accidental 
degenerades |3^,^]. A numerical plot of A(e) given in Fig.l shows that this term oscillates rapidly. The fitting of 
Sup e '< e A(e) gives Sup e '< e A(e) ~ e 7 where 7 = 0.6. Since 7 < 1, the first two terms in Eq. j5^ ) are still dominant as 
long as e >> 1. However, the contribution AÍx(e) arising from overcounting the points on coodinate axis is proportional 
to e 1 / 2 and smaller than the second term J5(J. Hence the fluctuating part of the cumulative density of states A(e) in 
Eq. (|5^) dominates over the contributions from A\ and Aq terms in Sec. 2.3. Similar arguments should hold in any 
finite size systems regardless of whether the system is integrable or not. 

For these reasons, here we properly take into account the lowest energy gap which carries essential information 
about finite size effects, and use the continuous spectrum approximation above the lowest excited mode. The density 
of states has the following form 

, \ n £l/2 n r 1 1 1 1 ,„> 

P 3 (e) = 7 T ïï + — + —] + •••• (53) 

4 aia 2 o aia 2 a\ a 2 

The first term is the Weyl term, the second term is the area contribution from the boundary. We can easily see these 
terms give the same terms in the heat kernel Eq. (^8|) related by a Laplace transformation. 
Next we write the heat kernel in terms of the density of states as 



K(t) = K + P3 (e)qlde, (54) 



3 Here we rewrite n as e (dimensionless energy). 
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where s\ = 1 corresponds to the energy level of the lowest excited mode and Kç, is the contribution in Eq. ( (49| ) from 
the ground state. Due to the presence of a cutoff, one can show that the total charge of all excited modes is 

Q = ò 2 T 2 + ^log^ + 0(A), (55) 
2 ui 

where m 2 = 7r 2 /L 3 2 + m 2 — fi 2 and O(A) is the contribution from the residual term A(e) and will be ignored hereafter. 
Then the second term in Eq. ( fÏ6| ) is replaced by b 5 / 2 T\ogTL^ for large L 3 close to the critical temperature. 

Now we evaluate the finite size correction to the critical temperature. The bulk critical temperature Tc is defined 

by 

Q = b 2 T^ 2 . (56) 



From Eq. (46) and the above argument, the leading correction to the bulk critical temperature for a finite system 
manifests as 

Q = b 2 T c 2 + b 3/2 T c logT c L 3 . (57) 
From Eq. ( |56] ) and Eq. (|57|), we obtain the finite number correction to the critical temperature for small b 3 / 2 

T c b 3/2 \og(QL 2 /b 2 ) 



T c (0) 4(6 2 Q)V2 

The correction to the condensation fraction can be easily obtained as 

Qo _ . _ Qi 
Q ~ Q 



T b 3/2 \og(QL 2 /b 2 ) T T 

&3/2 T T 

■ lo g(37nT)- 



(6 2 g)V2 Tc (o) *v T( (e». 

In Fig. 2, we plot the condensation fraction of the ground state as a function of the temperature. As mentioned in 
Section 2.3, the finite size correction in Eq. ( m ) gives the opposite shift in the critical temperature whcthcr Neumann 
or Dirichlet boundary condition is used. In the case of Dirichlet boundary condition, as we saw in Eq.(^3|), the surface 
term decreases the density of states. The smaller density of states requires the excitations with higher energy for 
the condensation criteria (|ï|) to be met and therefore the condensation has to start at a higher temperature. These 
results agree with those found in |ÏÏ],|5(J. It is of interest to compare the results obtained here with those for àtoms 
trapped in a harmònic oscillator potential. One can easily see that the boundary effect in such a potential is due to 
the Neumann boundary condition p0| ]: the surface term increases the density of states which results in the decrease 
of the critical temperature from the bulk value as observed in |33|. 



B. Multistep Condensation 



1. One-dimensional Condensation 



As mentioned in the Introduction, in the presence of strong anisotropy, condensation can occur in multisteps. To 
see one-dimensional condensation, we require L\ = L 2 << L 3 , equivalently, a± = a 2 >> 1 where a\L\ = a 2 L 2 = L 3 
as defined in Sec. 3.1. In such a case, it is meaningful to split all the excited quantum states into one, two, and 
three-dimensionally excited modes in the following way. Hereafter we focus on Neumann boundary conditions through 
the rest of the paper. 

The corresponding heat kernels for these states can be defined as 

00 

Ei 2 2 

n 3 = l 
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oo 

712 — 1 ri3— 

oo oo oo 

K 3 (t) = e- r 'i 2 " 2 "i 2r e-^ 2n * T e-^ n ' T , (60) 



ni— 1 n2 — 1 ri3=0 



respectively. The factor 2 in K 2 (t) is due to the symmetry between L x — L 3 plane and L 2 — L 3 plane. 

Following the similar steps from Eqs. ( J50| ) to (|53|), we obtain the expression for the three-dimensional density of 
statcs as 



P3(£) = T — - ^[72 +-]+•••■ ( 61 ) 



7T5 1 / 2 7T 1 2 

4 af 24 ^a 2 + ai 

And the three-dimensional heat kernel in terms of the density of states is given as 



/•OO 

K 3 (t) = / p 3 (e)q E 3 de, 

J S\ 

where e\ = a\. This gives the total charge of three-dimensionally excited modes 



(62) 



Qs = b 2 T 2 -^^log^, (63) 
6 m z 



where rh 2 = ir 2 Si/L 3 2 + m 2 — p 2 . 
For the two dimensional heat kernel, 



E2 2 2 % * 2 2 2 

e -r, 2 tt n 2 r ^ e -r, 3 tt n 3 r 

n 2 — 1 n 3 — 



00 00 

2^, 2 , „ 2 



2E E«3 in2+ " 3 = £2r 2 (n) g ™, (64) 

77-2 = 1 77-3—0 Tl— 



where ^(n) is the number of solutions of n = a^n^ + n 3 2 in natural numbers. We write Eq. (|64| ) in an integral form 
using the density of states p2^) — 7r/2ai such that 



K 2 {t) = / de P2 (e)q E 3 



ú 1 



2irair) 3 t 



(65) 



where £\ = a 2 for the present choice of units. Then from the expression of the zeta function in terms of the heat 
kernel in Eqs. ( p8| ) and (|2^) we obtain 

T T fÇ.s-2 rco 

Ca( S ) = ^W/ dr T s - 2 e-™ ^W^t/tt), (66) 



^ T(s) J 

where in 2 is the same as in Eq. (|6^) . The total charge of two-dimensionally excited states is 

2/iL 2 L 3 T T 2 

Q2 = log ■ (67) 

tt m z 

For the one dimensional case, we have 

52s /-oo 



ois r°o 2 

CaO») = £t / r- 1 2 gj» ^-^^(^rA). (68) 

1 W J ° «3=1 

The total charge carried by one-dimensionally excited states has the form 
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The three-dimensional critical temperature is reached when all the three-dimensionally excited modes are saturated, 
namely, sum of all the modes with energy larger than e = a\ is equal to the total charge Q. We write this condition 
corresponding to Eq. (0) as 

Q = Qa(T aD ). (70) 

Thus we obtain 

r\ urr2 h/2T 3D T 3D L X AmL 2 L 3 T 3D T 3D L X 

Q = b 2T 3D '— log + log , (71) 

ó 7T 7T 7T 

where we set [i = m. The third term is the contribution from two-dimensionally excited modes with energy larger 
than e = a 2 . We will ignore the contribution from Q\ based on the argument in Sec. 3.1 that this term is dominated 
by the residual term A(e). 

For one-dimensional condensation to be observable, we must have 

Tid < T 3D . (72) 



Furthcrmore for sufficiently large L 3 » Li,L 2 , comparing the expression for Qi in Eq. (69) with those for Q 2 and 
Q 3 in Eqs. (|67|) and (pa), we obtain 



Q2(T 1D ), Q 3 (T 1D ) « Qi(T 1D ). (73) 
Thus we obtain Tm as the temperature at which one-dimensionally excited states saturate, i.e. 

Q = Qx{Txd). (74) 

This gives T 1D = 7rQ[2rn J L 3 2 log(2 7 r)]- 1 . 



2. Two-dimensional Condensation 



For two-dimensional condensation, we assume L\ << L 2 = L 3 , whence we split the excited quantum states into 

oo 

Ki{t) = 2 ^2 e -n^W? 

n 3 — 1 

oo c 

E2 2 2 x ^ 
e ~ V2 * n 2 r J2 e 

n 2 — 1 n 3 — 1 

oo oo 

E 2 2 2 ^ ^ J J J ^ ^ J J j 



2 2 2 
-7J 3 7T U 3 T 



ii2=() n3— 



The integer a\ >> 1 defined by aiii = L 2 = L 3 will be used. The factor 2 in i^i(r) acounts for the symmetry 
between L 2 direction and L 3 direction. 

The three-dimensional density of states becomes 



n e 1 ' 2 tt_2 
4~af~ + 24 

This gives the total charge carried by three-dimensionally excited modes 



PÀs) = -— + 7T7[- + 1 ] + ···· (76) 



Q 3 = b 2 T 2 + h -llf\o g % (77) 



where rh 2 = ir 2 ei/L 2 + m 2 — fi 2 and s± = a 2 . 

The two dimensional heat kernel is given by the density of states p 2 {e) = 7r/4 as 
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K 2 (r) = dep 2 (e)ql = -4jL-, (78) 



where e\ = 1. The total charge of two-dimensionally excited states is 



pL 2 L 3 T T 2 

Q% = — = lo S^2- ? 9 



The total charge carried by one-dimensionally excited states is given by Qi = 4/iL 3 2 Tlog(27r)/7r. 

The three-dimensional critical temperature is reached under the same condition as in one-dimensional condensation, 
Q = Qs{T 3 d)- Thus we obtain 



2 b 3/2 T 3D T 3D Li mL 3 2 T 3D T 3D L X 



Q = b 2 T 3 f } + -^-^\og^^ + i-iüiog-ïí^, (80) 



where the third term is the contribution from the two-dimensional modes as in Eq. (|71[). 
For two-dimensional condensation to be observable, we must have 

T 2D < T 3D . (81) 

Furthermore, for sufficiently large L 2 , L 3 >> L%, by comparing the expression in Eqs. ( \Ï7\) and (|79|), we have 

Q 3 (T 2 d) « Q2(T 2D ). (82) 

Thus we obtain T 2 jj as the temperature at which two-dimensionally excited states saturate, 

Q = Q2(T 2D ). (83) 

Then T 2D = Q[b 3 / 2 \og(QL 3 /b 3 / 2 )}~ 1 for large Q, where b 3/2 = mL 2 L 3 /n = mL%/-n:. 

3. Three-step Condensation 

To show three-step condensation, we assume L\ << L 2 << L 3l or equivalently, a\ >> a 2 >> 1. The corresponding 
heat kernels for these states can be defined as 



Ki(r) 


OC 

v -s 2 2 2 

ns — l 


= Í[fl a (0|i»j3 2 7rr)-l], 






00 

X "V 22 2 


OO 1 


-l][e 3 (0|ÍJ73 2 7rr) + l], 




no — l 


n 3 =0 




K 3 {t) 


00 

v -s 2 2 2 

_ e -Vi 7T n x T 


00 00 
x — ^ 22 2 x — ^ 22 2 






7li=l 


n 2 =0 «3=0 






= \[e 3 (Q\iv^r) - 


-l][^(0|iry 2 2 7rr) + l][e 3 (0|í% 2 7rr)- 





(84) 



The asymptotic behavior when 1 >> r\\ >> 772 >> V3 can be derived in a similar way as in Eq. ( |38| 
The three-dimensional density of states becomes 



/ n T eV2 T r 1 1 1 1 

Ps(e) = j + 7T A \ + — + —] + ••• (85) 

and the total charge of three-dimensionally excited modes is given by the same form as in Eq. (^) . 

From the two dimensional heat kernel, we obtain the two-dimensional density of states p 2 (s) = 7r/4a2 such that 



K 2 (t) = / de P2 (e)q £ 3 

Jsi 



4ira 2 r) 3 r 
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where E\ = a 2 for the present case. This will give us the total charge of two-dimensionally excited states as 



Q 2 = ^^log^. (87) 



liL 2 L 3 T T 2 
-7. log — 

The total charge carried by one-dimensionally excited states has the same form as in Eq. ( p9| 
The three-dimensional critical temperature is obtained as 

2 . fa/iTsD T 3D Li mL 2 L 3 T 3D T 3D Li 



Q = b 2 T 3 z D + -^-^fog-^ + ^^log 

à 7T 7T /l 

To observe three-step condensation, we must have 

T 1D < T 2D < T 3D . (89) 

In addition. for large anisotropy L\ << L 2 << L 3 , comparison of explicit formulas for Qi,Q 2 and Q 3 in Eqs. (|6S|). 
© and (0) gives 

Q 3 {T 2D ) <<Q 2 (T 2D ) (90) 

and 

Q 3 (T 1D ) « Q 2 {T 1D ) « Qi(Tid), (91) 

where T 2 d (T\d) is obtained by the saturation of the two (one)-dimensionally excited states as Q = Q 2 {T 2 £>) and 
Q = Qi(Tid). T 2 d is given in the leading order by 

Q = b 3/2 T 2D \og{L 2 T 2D ) 1 (92) 

where b 3 / 2 = mL 2 L 3 /ir. We thus obtain T 2 o = Q[b 3 / 2 log(Q L 2 / 'b 3 / 2 )]~ l for large Q. T\u has the same form as in the 
one-dimensional condensation case. 

From Eqs. (|9|) to (|9Ï|), we obtain the following inequalities for anisotropy parameters and the total charge corre- 
sponding to (A) T 1D < T 2D , (B) T 2D < T 3D , and (C) T 1D < T 3D as 

h » lQ g^) (93) 
L 2 21og(27r)' [ ^ } 

L3 » À> (94) 



Li 3[logQ] 2 



L·l 
Li 



L 



» i^ïïoi^F' ^ 



where Q = TrQ/mL 2 . In Fig. 3, different multistep behaviors corresponding to various ranges of anisotropy parameters 
are shown. Three-step BEC can be seen in a wide region of parameter space. For an extremely strong anisotropy, 
dynamics along the short edge length will freeze out (EIRD < 3) before BEC into Q 2 or Qi sets in. In such a case, 
we only observe two or one-step BEC. For a quasi-linear cavity along the vertical axis, three-step BEC can still be 
observed, while for a quasi-planar cavity along the horizontal axes, only up to two-step BEC can happen. 

In Fig. 4, the condensation fractions Qo/Q,Qi/Q,Q 2 /Q,Q 3 /Q as a function of the temperature are plotted. In 
the isotropic case (Fig. 4a), condensation is only into the ground state. Due to the finite size effects, condensation 
occurs before the critical temperature is reached. In strongly anisotropic cases, condensation occurs in steps. In Fig. 
4b, one-dimensional condensation is seen. T 3 d determines the onset of condensation into one-dimensionally excited 
states. Note that at T 3 d, the ground state fraction is negligibly small. Condensation into the ground state occurs at a 
much lower temperature. In Fig. 4c, two-dimensional condensation is plotted. At T 3 d, two-dimensional condensation 
manifests itself. The critical condition fi = m is satisfied well in Fig. 4a-c. In Fig. 4d, three-step condensation is 
shown. Three-dimensionally excited modes dominant in higher temperature are condensed into two, one, and the 
ground state as the temperature is lowered. The deviation of T 3 o and the onset of two-dimensional components 
reflects the fact that the condition fj, = m is not satisfied for the parameters chosen in Fig. 4d at T 3 o- This is another 
manifestation of finite size effects and the result should improve near the thermodynamic limit (rji —> 0, Q — > oo). The 
similarity between each condensation process becomes evident in the logarithmic T scale as can be seen in Fig. 4e. 
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In conclusion, finite size effects on the Bose-Einstein condensation of a charged scalar field can lead to the multistep 
condensation in the presence of strong anisotropy. 

In this paper, we started from calculating the effective action to one-loop order using zeta function regularization. 
Large volume and small mass conditions are assumed to facilitate an asymptotic expansion of the heat kernel and 
finite size corrections corresponding to the surface term, còrner term, etc. are obtained. We proceeded beyond the 
continuum spectrum approximation and showed that the higher order terms in the Standard asymptotic expansion 
are dominated by the contribution from the fiuctuating part of the density of states due to accidental degeneracy. 
The lowest energy gap is shown to play the crucial role in determining the critical temperatures for one and two- 
dimensional systems. The corresponding low-dimensional critical temperatures are calculated. The energy spectrum 
and the associated heat kernel can be partitioned into parcels of eigenmodes excitable in dimensions 3, 2, 1, or 0. As 
the temperature is lowered, modes in different parcels behave quite differently in the presence of strong anisotropy. 
Whcn Ti£,,T-2D < T3D are satisfied, condensation occurs first into the lower dimensionally excited states at T^d 
following the ground state condensation at lower- dimcnsional critical temperature. Experimental observation of 
these phenomena can in principle be realized in an anisotropic harmònic potential traps. 
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In this appendix, the exact formula for Aíd(s) for arbitrary dimension d is derived. <i-dimensional cumulative 
density of states A/d(e) can be obtained straightforwardly from Àíd(s) as we showed in Section 3. Suppose Aíd{£) 
counts the number of integer solutions n of the equation e = a 2 n 2 + ■ • ■ + a d 2 n d 2 where a = (01, • • • , ad) is a constant 
cí-dimensional vector with integer coordinates. Then J\fd(s) can be written as 



APPENDIX A: CUMULATIVE DENSITY OF STATES 






(Al) 



where A = a% ■ ■ ■ ad, an = (dini, • • • , a^n^), and l/a = (Zi/di, • • • , ld/ad)- Summation is over all <i-dimensional vectors 
with interger coordinates. Poisson's summation formula is used to obtain the third line and 




(A2) 




This yields 




for 1^0 
for í=0, 
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where Jd(x) is the Bessel function and Va = ir d / 2 /T{d/2 + 1) is the volume of a d-dimensional sphere with unit radius. 
Hcnce we obtain 

M d{ e) = A-^V^ + A-W £ Jd ^^f ) - (A4) 

\l/a\ a >' s 
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Figure Captions 



Figure 1 The residual term A(e) in the cumulative density of states (Eq. (]5^)) is plotted for the Neumann boundary 
condition. This term shows the highly oscillating behavior due to accidental degenerades. Fig.la is an isotropic case 
(ai =«2 = 1) and Fig.lb is an anisotropic case (a\ = 10 and et2 = 3). Supremum (dashcd lines) of both curves show 
that the rate of increase is proportional to e 7 where 7 = 0.6. 



Figure 2 The condensation fraction Qo/Q for the Neumann boundary condition is plotted as a function of the 
temperature. L\ = 1, L% = 10, L3 = 100, Q = 10000, and m = 0.1. Dotted curve shows the bulk contribution. 
Sòlid curve includes the finite size correction based on Eq. (|9|) . T c denotes the finite size corrected three-dimensional 
critical temperature defined in Eq. (S). 

Figure 3 Different multistep behaviors corresponding to different anisotropy parameters L2/L1 and L3/L2 are indi- 
catcd. The logarithmic scale is used for both axis. Q = nQ/mL2 — 10000 is fixed. Multistep BEC can only happen 
in an intermediate yet strongly anisotropic regime. 



Figure 4 The condensation fractions Qo/Q (sòlid curve), Qi/Q (dashed curve), Q2/Q (dot-dashed curve), Q3/Q 
(dotted curve) as a function of the temperature are plotted for the Neumann boundary condition. Isotropic case 
(Li = L2 = L% = 3, Q = 100, and m = 2) are plotted in Fig. 4a. Condensation is only into the ground state. 
T c = 1.97 is the critical temperature in Eq. (|57|). In Fig. 4b-e, anisotropic cases are shown. In Fig. 4b, L\ = 2, 
L2 = 2, L3 — 300, Q = 2000, and m = 1 are chosen. One-dimensional condensation occurs in this case. T c = 2.03 
is the three-dimensional critical temperature in Eq. (|7l|). In Fig. 4c, L\ — 2, L 2 = 200, L% = 200, Q = 8000, 
and m = 0.5 are chosen. T c = 0.98 is the three-dimensional critical temperature in Eq. (|80|). Two-dimensional 
condensation can be seen. In Fig. 4d, condensation occurs in three-steps. L\ = 2, L 2 = 100, L 3 = 600, Q — 4000, 
and m = 0.5 are used. The long dashcd line is the chemical potential /i. T c = 0.79 is the three-dimensional critical 



temperature in Eq. (88). The logarithmic T scale is used in Fig. 4e for the parameters in Fig. 4d. 
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